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ABSTRACT 

Photographs of 521 shower cores in an array of current-limited spark 
( 'discharge 1 ) chambers at Sacramento Peak' (2900m above sea level, 730 g 
cm“ 2 ) , New Mexico, U.S.A., have been analysed and the results compared 
with similar data from Leeds (80m above sea level, 1020 g cm” 2 ) . It was 
found that the 'central' density differential spectrum is consistent 
with a power law index of -2 up to ^1500nT 2 where it steepens, and that 
shower cores become flatter on average with increasing size. 

Scaling model predictions for proton primaries with a «e”2.71 ener gy 
spectrum account well for the altitude dependence of the data at lower 
densities. However, deviations at higher densities indicate a change 
in hadron interaction characteristics between ^ f ew x 10 14 and 10 15 eV 
primary energy causing particles close to the shower axis to be spread 
further out. 


1. Introduction . Experimental data on particle distributions within a 
few metres of the axes of air showers at Sacramento Peak (2900m above 
sea level, 730 g cm -2 ), New Mexico, U.S.A., are described. The data are 
the results of the author’s analysis of photographs of showers in an 
array of current-limited spark chambers (’discharge’ chambers) operated 
by Hazen et al. (1981). These photographs were originally obtained as 
part of a search for high-pj- subcores. 

The results are compared with similar data from the 35m 2 discharge 
chamber array at Leeds (80m above sea level, 1020 g cm -2 ) (Hodson et al., 
1985) and predictions from Monte Carlo simulations. 

2. Experimental Details . The array was housed under a canvas tent and 
was photographed from above using a small camera. During the experimen— 
tal runs used as the basis for the present work the array consisted of 

20 lm discharge chambers (nearly identical to those used at Leeds) 
arranged into A parallel rows each of 5 chambers, with -v30cm gaps (walk- 
ways) between adjacent rows and the chambers in each row close-packed to 
give in most cases only ^3cm width of dead space between them. The 
trigger condition was a particle density *50nT 2 in a scintillator near 
the centre of the array in coincidence with at least one particle in a 
scintillator a few metres outside the array. The trigger rate was 
^I0hr“ 1 . 

3 « Analysis and Results. The methods of analysis of the shower 

photographs follow closely those of Hodson et al. (1983a, b) for the Leeds 
data. The photographs (negatives on 35mm film) were projected at a 
linear demagnification of 8 in array space and scanned by eye for cores. 
Core location and subsequent density measurements made use of an annular 
grid, superposed over the projected image, identical in pattern to that 
used for the Leeds photographs (see Hodson et al., 1983a) but with its 
dimensions (in array space) linearly scaled up by a factor s=(1020/730) . 
This was to take approximately into account the effect of the difference 
in air density on particle densities when comparing the data for Leeds 
and Sacramento Peak. As in Hodson et al. (1983a) the grid was used to 
determine a shower centre (centre of symmetry) for the particle distri- 
bution of each core located in the array area and, with the grid 
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centred on the shower centre, densities at various radial distances were 
determined. To reduce the effects of trigger bias and scanning bias the 
analysis excluded all showers with centres located in the walk-ways 
between the rows of chambers, those with centres falling in the 4 corner 
chambers and any in the remaining outer chambers with centres falling 
within 0,4m of the edge of the array. The total collecting area for 
shower cores was then equal to -12m . 

A count of discharges in the centre circle (radius 0-25sm) of the 
grid gave the f central 1 density; when a walk-way overlapped part of 
the centre circle (maximum 50%) then symmetry was used to obtain the 
count for the whole of its area. Counts in the angular sections of the 
annuli at the scaled distances sm and 2*5sm gave densities at these 
distances. (In general, counts were made in only one randomly-selected 
section in each of these annuli.) 

In the rest of this paper, observed particle densities per m 
(i.e, count/actual bin area in array space) are given, rather than 
densities scaled up by a factor s 2 . 




(a) (.b) 

Figure 1 Differential density spectra p(r) observed at Leeds (L) and 

Sacramento Peak (SP) 

Figure 1 shows the differential density spectra for 521 cores at 0, 
lrOs, and 2*5 s m from the shower centre (normalised to a 35ra 2 collecting 
area) as observed at Sacramento Peak, compared with the corresponding 
experimental data for Leeds, Power law lines of index -2 have been super- 
posed on the spectra to show the steepening at high densities in the Leeds 
data (see Hodson et al., 1985) and for comparison with Sacramento Peak. 
(The deviations of the spectra from the -2 power law at low densities are 
due to trigger and scanning bias.) In Figure 1(a) there is evidence of a 
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Figure 1(c) 


steepening from the -2 power law in 
the Sacramento Peak p(0) spectrum at 
' v 'l * 5 . 10 3 m -2 ; the spectra further 

out (Figures 1(b) and 1(c)) may also 
steepen around 10 3 m~ 2 but this is 
less clear. 

Figure 3 gives a plot of core 
'flatness', <p ( 1*0 s) >/<p (0) > , versus 
P (2*5 s) (used as a measure of shower 
size, as p(4*0) in Hodson et al., 1985); 
this shows a tendency for shower cores 
to become flatter within s(=l # 4)m as 
the showers become larger (as also 
reported for the Leeds data (Hodson 
et al., 1985)). Taking Figures 1 and 
3 together, steep dense shower cores 
are observed at a lower rate than 
might be expected by extrapolation. 

Figures 2 and 3 give comparisons 
with simulation predictions, after 
excluding all real and simulated 
showers with p(0)<100nf 2 , using scaling 
models (with EMC Model I) for proton 




(a) # (b) 

Figure 2 Comparison of observed differential density spectra 
(histograms) at Leeds (L) and Sacramento Peak (SP) with simulation 
predictions (points) 

primaries, described oy Ash (1985a). The full detailed core-mapping 
calculations were used for comparison with the Leeds data; for 
Sacramento Peak a simpler calculation was used, with the shower centre 
assumed to coincide with the shower axis. For the present purposes this 
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Figure 3 'Core flattening' 
at Sacramento Peak 

difference is unimpqrtant. The proton 
primary energy spectrum was 
dn/dE = 1*36 x 10 4 E~ 2 ^* 71 Gev^ 1 m -1 sr -1 s" 1 , 
consistent with balloon-based observa- 
tions, here extrapolated to 10 16 eV. 

The predictions agree well with the observed spectra both in slope and 
in absolute rates in the -2 power law region, where proton showers would 
be expected to dominate the spectra. At higher densities, all the Leeds 
spectra and the Sacramento Peak p(0) spectrum are clearly steeper than 
those predicted. * 


Figure 3 shows that the predicted large proton showers have much 
steeper cores, on average, than those observed, as was also found for 
the Leeds data (Ash, 1985b). 

4. Discussion * It is interesting to consider the mean primary energy 

<E> of showers with a given p(r) at Sacramento Peak. At p(0) - 1*5 * 10 ^ 

m** 2 the proton simulations predict <E> -6x 10 14 eV; if the steepening of 

the p(0) spectrum is due to a rapid increase in the slope of the proton 
primary energy spectrum then a clear steepening should also occur for 

similar values of <E> in the density spectra further out. At 2*5s and s m 

from the shower centre <E> = 6 x 10 14 eV corresponds to densities of 'V' 200m 2 
and 400m” 2 respectively; no steepening is evident at around these 
densities. This, together with the observed core flattening, lends 
support for a significant change in hadron interaction characteristics 
between primary energies of ^ a f ew x 10 14 and 10 15 eV. A natural expla- 
nation might be found in a mechanism which spreads particles out from 
within ^lm of the shower centre. Complete understanding probably 
requires the combined effect of interaction changes plus some steepening 
in the proton energy spectrum. 
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